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Abstract —Diffusion-based and nenral communication are two 
interesting domains in molecular communication. Both of them 
have distinct advantages and are exploited separately in many 
works. However, in some cases, neural and diffnslon-based ways 
have to work together for a commnnication. Therefore, in this 
paper, we propose a hybrid commnnication system, in which 
the diffnsion-based and neural communication channels are 
contained. Multiple connection nano-devices (CND) are nsed to 
connect the two channels. We define the practice function of the 
CNDs and develop the mechanism of exchanging information 
from diffnsion-based to neural channel, based on the biological 
characteristics of the both channels. In addition, we establish a 
brief mathematical model to present the complete communication 
process of the hybrid system. The information exchange process 
at the CNDs is shown in the simulation. The bit error rate (BER) 
indicator is nsed to verify the reliability of communication. The 
resnlt reveals that based on the biological channels, optimizing 
some parameters of nano-devices could improve the reliability 
performance. 

I. Introduction 

Molecular communication (MC), one of the most promis¬ 
ing communication ways at nanoscale, develops rapidly due 
to great advance of the nanotechnology. The information is 
encoded in chemical molecules, transmitted via biological 
channels and exchanged between bio-inspired nano-devices 
[1]. In general, molecular communication could be divided 
into wireless and wired branches. Wireless MC is realized 
through diffusive propagation of the molecules, and could be 
further divided according to the communication distance and 
molecule species. Some typical examples includes bacterium, 
calcium signaling and pheromone. In wired MC, information 
is transmitted along some types of physical link, such as the 
neuron, blood vessel and microtubule. 

Diffusion-based communication is an important branch 
of wireless MC, in which the molecules could only trans¬ 
fer slowly in a short distance. While nano-devices adopting 
diffusion-based communication could move freely in the envi¬ 
ronment rather than stay at a fixed position. Differently, neural 
communication, one typical example of the wired MC, is a fast 
and long range transmission way in human body. However, 
nano-devices of neural system lack flexibility, i.e., they have 
to stay at some key position to proceed the communication. Eor 
example, receiver should be close to the cell membrane, for 
monitoring the membrane potential during the communication. 
In nature, there are some cases that diffusion-based and neural 
ways cooperate to proceed the signal transmission, and we list 
two cases. Eor the first case, the Ca^+ ions diffuse around 


the neurons, access the cell through Ca^'*' channels, impact 
the vesicle release activity to adjust the neural signaling. For 
the second case, some ions including Na+, K+, Cl“, diffuse 
and pass in and out the neuron through ion channels, altering 
the membrane potential to trigger the electric signal of neural 
fiber. In this paper, to share the different advantages and avoid 
shortages of the both communication ways referred before, 
we design an adjustable diffusion-neuron hybrid system by 
extending the biological background of the second case. 

From the aspect of communication and mathematical 
scheme, abundant literatures investigate the two communi¬ 
cation ways discretely, which lay a theoretical basis for our 
research. Such as [2], [3], [4], respectively discuss the channel 
capacity, transmission reliability and delay for controllable 
diffusion-based MC. About the neural communication, in [5], 
a mathematical model of single neuron at Cornu Ammonis 
region is established. [6] analyzes the error probability and 
delay of neural transmission. In [7], transmission interference 
in neural channel is studied. 

However, with the aid of those basis, it’s still difficult to 
establish a hybrid communication architecture. In traditional 
communication, a gateway device is necessary for connecting 
the wired and wireless networks. An typical example is the 
home router, with wired and wireless interfaces integrated in 
one device, coping with the information exchange between 
different mediums. At nanoscale, a similar device is needed 
to realize the same function in the hybrid architecture, which 
is named as connecting nano-device (CND) in this paper. The 
major challenge is a reasonable design of the CND, and a 
feasible mechanism to exchange signal from diffusion-based 
channel to neural channel. 

In this paper, we propose a diffusion-neuron hybrid com¬ 
munication architecture, in which the information is transmit¬ 
ted from diffusion-based to neural channel, not including the 
reverse process. Different from [5]-[7], in our work, we con¬ 
sider the general controllable diffusion process between nano¬ 
devices, not limited by the ion diffusion of the natural neural 
activity. Due to the restriction of size and power for nanoscale 
device, a CND cluster is utilized rather than a single CND to 
perform the task. We define two interfaces for CNDs, named 
diffusion interface and neural interface, respectively coping 
with the communication for different mediums. We design the 
practice function of the CNDs, which apply the electric current 
stimuli to trigger the neural signaling. This method is utilized 
in [8], [18], [19] and [20], for design of the nano-device in neu¬ 
ral communication. [8] proposes the equivalent neural devices 
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Fig. 1. Biological model of the hybrid communication system. 


method, that the multiple neural devices are uniformly driven 
by one current. While in our paper, for the neural interface of 
each CND, we apply the different strength of currents, that is 
inspired according to signaling of the diffusion interface. We 
describe it as the current inspired mechanism, through which 
the signal is exchanged from diffusion-based to neural channel. 
We discuss this mechanism in biological method and present 
it in the complete communication process with mathematical 
model. The expression of bit error rate (BER) is given, to 
study the reliability of the hybrid model. In the simulation, 
we study how the information is exchanged at the CNDs, as 
well as which parameters impact on the BER performance. 
Our work provides some guidance for designing the hybrid 
communication architectures at nanoscale, contributing to the 
development of body area network. 

The rest of the paper is organized as follows. In Section 
II, we introduce the biological background of diffusion and 
neural communication. Based on this, the entire hybrid system 
is described, the CND scheme is discussed. The mathematical 
model of a complete communication process is established in 
Section III. After that, the simulation result is presented in 
section IV. Einally, we conclude the paper in section V, as 
well as the future extension of this work. 

II. Biological Method 

Natural diffusion phenomena exist abundantly around the 
neural system. Such as the ions Ca^'*', Na+, K'*", Cl“, diffuse 
around the neurons, and play an important role in normal 
activity. Though they are parts of the neural signaling, the 
diffusion behaviour could be extended into general diffusion- 
based communication by applying more types of molecules 
in artificial nano-devices, to implement the interface of neural 
network to other MC networks. The neural transmission is 
rather complex, proceeded among billions of neurons. The 
typical communication process in a single neuron includes 
neural firing, action potential (AP) transmission, vesicle release 
and postsynaptic response. Briefly, information is encoded in 
AP spike trains, which are generated in process of neural firing. 
The potential rises and falls rapidly like water wave, moves 
along the neural fiber until aiTiving presynapse, as Pig. 1 (b) 
shown. Then, ion channels are open, letting Ca^+ ions come in 
and promoting the process of vesicle release. Neurotransmitters 
stored in vesicles are released, and propagate in the synaptic 
gap before arriving in postsynapse, where the postsynaptic 


response happens to raise the postsynaptic membrane potential, 
due to the exchange of various ions. Postsynaptic response po¬ 
tential is divided into excitatory postsynaptic Potential (EPSP) 
and inhibitory postsynaptic Potential (IPSP), both of which are 
much lower than one AP. The former contributes the neural 
firing process while the latter restrains. With the effect of 
EPSPs and IPSPs, when membrane total potential exceeds a 
threshold in postsynaptic response, an AP spike is generated, 
followed by a period of resting time, named absolute refractory 
period in neuroscience. During the absolute refractory period, 
no AP spikes could be generated. 

The hybrid system designed is composed of one transmitter 
(TX), one receiver (RX) and multiple connecting nano-devices 
(CND), across both the diffusion-based and neural communi¬ 
cation mediums. The biological model is shown in Pig. 1 (a). 
TX is a biological nano-device as most literatures desire, which 
could produce, emit and receive specific type of molecules. RX 
considered in this paper is the typical nano-device in neuron 
communication, with ability of monitoring membrane potential 
and decoding information. CND cluster are designed as the 
relaying bridge, which could receive molecules from diffusion- 
based channel and evoke EPSPs to fire the AP spike trains in 
the neural channel. 

As a key element of the hybrid system, CND cluster should 
be able to cope with the communication both in diffusion- 
based and neural channel. The CND cluster around a neuron 
is shown in Pig. 1 (c). They are bio-inspired nano-devices 
containing diffusion and neural interface. Responsible for 
communicating with TX, the diffusion interface is formed 
by plenty of biological receptors, which are embedded on 
the membrane surface, distributed uniformly to receive signal 
molecules in a small receiving space 14- We utilize the ligand 
receptors for diffusion-based interference described in [9]. 
Molecules are bound and released with different rates, when 
they are in the receiving space. We adopt the design of [8] 
for the neural interface. CNDs should be able to produce tiny 
biological electricity to stimulate the neural firing process. 
An example of generating biological electricity is introduced 
in [17], depending on oxidation of organic substrates, which 
is promising to implement in an artificial cell. The neural 
membrane and CND are directly contacted with gap junctions, 
which is also the connection way of two adjacent neurons. 
Clusters of connexine proteins (named Cx36), combine to form 
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Fig. 2. Communication process diagram of the hybrid system. 


an ion channel, allowing the ions and tiny electricity passing 
through during the neural encoding process [10]. The gap 
junction method is supported by neuroscience in [11] and [12], 
where it is studied in oscillatory behaviors and synchronization 
phenomena between neurons. 

The on-off key (OOK) modulation of diffusion-based MC 
is used at TX, determining that CNDs could receive much 
molecules for bit ”1” and few molecules for bit ”0”. The idea 
of current inspired mechanism is that according to quantity 
of the received molecules per bit, CNDs operate indepen¬ 
dently and inspire the electric currents of various strengths, 
to implement the device-based neural firing process. First of 
all, to avoid the possible interference, this process should 
proceed when there is no normal neural activity. Result of 
[18] indicates that different shape of current waves make 
difference for neural firing and the sine wave is the best. 
The current of several p.Ajcmr'^ is enough to fire an AR 
Hence, in this paper, we consider the sine wave current with 
various amplitudes. During this process, there is no decoding 
or complex calculation for CNDs. The sine wave amplitude 
just varies with the received molecules quantity per bit in 
different time slots. Then, different strength of currents evoke 
the EPSPs with various amplitudes, which contribute to raise 
the membrane total potential. As described before, the AP 
spikes could be generated only when membrane total potential 
is strong enough to exceed the threshold. Hence, the binary 
sequence information is encoded in the AP spike trains (no 
spikes = ”0”, spikes = ”1”), transmitted in the neural network 
until decoded by RX. 

HI. Mathematical model 

To detail the design of the hybrid system described above, 
we present the communication process from TX to RX with 
a brief mathematical model. As Fig. 2 shown, the communi¬ 
cation process is divided into 5 parts, which are transmitter 
encoding, diffusion process, CND operation, neural transmis¬ 
sion and receiver decoding respectively. We describe it as 
follows: an input binary information sequence, denoted by 
XiA = 1,2,..., is encoded in molecules by TX and emitted 
into the medium. Those molecules diffuse freely until received 
by CNDs. Then, CNDs stimulate the membrane of the neuron 
cell through tiny electronic currents, firing the AP spike trains 
S{t), in which relaying binary sequence Yi is coded. After a 
distance of neural transmission, the AP trains are decoded to 
Zi by RX through monitoring the membrane potential. 

A. Transmitter Encoding 

In this paper, we only consider one transmitter. One time 
slot is utilized to transmit one bit, with fixed length of T. 
Let R{0,t) be the molecular emitting rate of transmitter, at 


location of 0 and time of t. The OOK modulation is adopted 
to encode the input signal Xi, meaning that at the beginning 
of each time slot, quantity of Qtx molecules are emitted to 
indicate bit ”1”, i.e.. 


R{0,t) 


fQtx bit’T” 
\0 bit ”0” 


( 1 ) 


B. Diffusion Process 

Molecules diffuse in the environment and perform the 
brownian movement. For the diffusion process of single source, 
the molecular concentration follows the Fick’s second law, 
with the solution of the Green function [1], 

c{x,t) = ■ g{.xA)+ni{t) (2) 

Where c{x, t) is the molecular concentration at location x and 
time t. ni{t) is the gaussian noise of diffusion channel, s.t., 
ni{t) ~ N{0,a1). g{x,t) is the Green function mentioned 
above, expressed as g{x,t) = ea;p(-j^). Here D 

is the diffusion coefficient related with molecule species and 
environment condition. When t = kT, there have been k bits 
sent from the transmitter. Let r be the distance between TX 
and a CND, the concentration at a CND is accumulated with 
time and expressed as, 

k 

c{r, t) = f?(0, t — iT) ■ g{r,t — iT) + ni (f) (3) 

i=l 

C. CDN Operation 

Assuming there are M CNDs in the cluster. Let r*, be 
the distance between transmitter and CND. Utilizing the 
ligand-based receptors for diffusion interface, according to (3), 
for bit at CND, we calculate the quantity of received 
molecules as. 


- V ■ D 

Qrxik,i) = - - -/ c{rk,t)dt (4) 

e-l JiT 

Where ei and e_i respectively means the binding and release 
rate of the molecules [9]. p is the concentration of the ligand- 
based receptors on surface of CNDs. Vr is the receiving space 
of one CND. Based on the current inspired mechanism in 
section II, electric current Af-it), inspired by k*^ CND, are 
produced according to knowledge of Qrxik,i)- So, we have, 

Ak{t) = fiQrxik, t)), iT<t<{t + 1)T (5) 

Where f{x) is a positive relation function of x, indicating the 
current inspired mechanism. During the time slot, more 
molecules received, stronger the electricity should be inspired. 
Applying sine wave form, we rewrite (5) as, 

Ak{t) = p ■ Qrx{k,i) ■ sin{at Yb) (6) 

Where p indicates the sensitivity of the current inspired 
mechanism, meaning the current strength per molecule, a is the 
fundamental frequency of the sine wave, b means the phase, 
which is usually set as 0. 

With the given electric currents from all CNDs, the post- 
synaptic potential V(t) could be obtained according to Spike 





















Fig. 3. Transmission diagram of a single neuron 


Response Model, described in [13], 


E. Receiver Decoding 

The decoding process is monitoring the membrane poten¬ 
tial to obtain the AP spike trains S{t). On general condition, 
neural transmission is an efficient way with extremely low 
error probability because neurons has ability to eliminate 
interference [15]. So, S{t) could be seen undistorted when 
decoded by RX. Described in [16], it is efficient to compare 
a threshold 02 to correlation 7 , between membrane potential 
and h[t). The larger probability to inspire the APs, the bigger 
7 is. For time slot, iT <t<(i + 1)T, output Zi could be 
determined as. 


AT rtoo 

y{t) = yo + '^ i^{s)Ak{t - s)ds (7) 

k 

Where Vq is the resting potential without any neural activity. 
The current of CND is generated at time of tk- v is the 
linear response of the cellular membrane to an input pulse of 
current, given that the neuron did not hre in the recent time. 
Akit) means the current fired by CND. The AP spike trains 
are generated only when V(t) is strong enough to exceed the 
firing threshold 0i, the typical form of the AP spike trains is 
given by. 


f 1 7 > 02 
7 < 02 


( 10 ) 


The reliability of the system is verihed by indicator of 
the bit error rate (BER). Based on the communication process 
diagram in Fig. 2, we give the BER expression of the hybrid 
system as. 


Pe=pr{Z, = 0\X, = l)+pr{Z, = 1\X, = 0) (11) 

Considering the middle binary sequence Yi, we have. 


sh) = ( 8 ) 

j 

Where — tj) means the shape of one spike, which is 
triggered at time tj if V{t) = 0i and > 0. It reaches 

the maximum after a short time and quickly drop, followed by 
absolute refractory period of over 15 ms, preparing generating 
the next AP spike. During time slot, if at least one AP 
spike ip exists in S{t), Yi is ”1”, else Yi is ”0”. 


D. Neural Transmission 

Inspired by CNDs, signal should pass plenty of neurons 
before decoded by RX. The transmission could be regarded 
as the repeated process of one single neuron transmission. A 
mathematical model of a single neuron communication system 
is established in [5], shown in Fig. 3. m{t) is the stimuli 
signal for neural firing, regarded as the input. LNP model 
contains the process of neural firing and AP transmission, that 
is introduced in section II and illustrated detailedly in [14]. 
S{t) is the AP spike trains transmitted in this neuron. Vesicle 
release is a probabilistic model. Neurotransmitters release from 
the vesicles, with probability of p, which is related with the 
neuron types. In postsynaptic response block, the random value 
q, with the probability density P{q), indicates the variability 
of the EPSP amplitudes. EPSP shape is the second step of the 
postsynaptic response, moulding the voltage by alpha function 
h{t), expressed as, 

h{t) = ^ ■ t ■ exp{l - (9) 

Where hp is the EPSP magnitude and tp indicates correspond¬ 
ing time of magnitude peak. n 2 {t) is the neural noise in the 
postsynaptic potential, modeled as additive gaussian noise, s.t., 
n 2 {t) ~ A^(0,cr|). The output signal z{t), which is also the 
stimuli signal of the next neuron, is formed at the postsynapse. 


Pe =pr{Y, = 0|X, = l)pr{Z, = 0|F, = 0) 

+ pr{Y, = 1|X, = l)pr{Z, = 0|y, = 1) 

+ pr{Y, = 0\X, = 0)pr{Z, = l\Y, = 0) ^ ’ 

+ pr{Y, = 1\X, = 0)pr(Z, = 0|y, = 1) 

Note that pr{Yi\Xi) is related with TX encoding, diffusion 
process and CND operation. Differently, pr{Zi\Yi) is related 
with neural transmission and receiver decoding. We could see 
that for the diffusion and neuron sides, if information of the 
two sides are transmitted both accurately or both wrong, then 
the final decoding of receiver is right. Else, the receiver will 
get the error bit. 

IV. Numerical Result 

In this section, we simulate diffusion-neuron hybrid sys¬ 
tem detailed in section II and III. We particularly focus on 
the signal transform process from diffusion-based to Neural 
channel. Some controllable nano-device setting parameters, 
including CND quantity M, CND sensitivity p, CND funda¬ 
mental frequency a, time slot T and receiver threshold 02 , are 
investigated in terms of their impact to BER performance. The 
goal of the simulation is to test the signal exchange process 
in CND operation, and show that optimizing those parameters 
could improve the communication reliability. 

A. Parameter Setting 

The results are computed under a common set of parame¬ 
ters, which are set as follows. The diffusion-based coefficient 
D is set as QApm'^/s. Molecule quantity of bit ”1” Qtx is 
fixed as 10^ [3]. Length of one time slot T is set from 100 to 
300ms. Distances between TX and different CNDs is set in 
range of 15-20 pm. ei and e_i are set as 0.1 and 0.08. p is set 
as 0.5 /xmol/liter, similar with [9]. The quantity of CNDs M 
changes from 5 to 12. electricity per receiving molecule p is 
set from 20 to 40 nAlcm~‘^ per molecule. The fundamental 
frequency a changes from 40 to 80 Hz. b is fixed as 0. The 
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Fig. 4. Neural firing for various bits, inspired by multiple CNDs, 
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Fig. 5. (a) molecular concentration c(rt^,t), (b) electric current Ak(t), (c) 

AP spike trains S(i), for bit ”1011”, T = 100ms, = 15^m, 02 = 1,M = 

6, )jL=20nAIcm~'^ per molecule. 


initial resting membrane potential Vo is set as -70mV for a 
typical neuron, with corresponding firing threshold 9i as - 
55mV. Receiver decoding threshold 02 is set as 1 similar with 
[4], EPSP magnitude hp and time to peak tp is respectively 
ImV and 0.5ms when calculating correlation 7 . The variance 
of diffusive noise ui and neural noise CT 2 are set as 0 . 1 //rm^ 
(average one molecule in 10 /im^) and QAmV. 

B. CND Operation 

We simulate the CND operations in computer according 
to the well known Hodgkin Huxley model [13], and present 
the signal exchange process from the diffusion-based to neural 
channel, which is the key of the hybrid system. 

In Fig. 4, We show the neural firing process at the CDNs. 
Each small wave curve means one EPSP contribution from 
CNDs. They have various amplitudes, decided by quantity of 
the received molecules. We can see that the potentials vary 
from -70mV, for bit ”1” and ”0”, there is an obvious difference. 
Membrane Voltage of bit ”0” increase slowly in Fig. 4 (a), the 
increase speed even falls behind the EPSP decreasing. So the 
potential could not reach 6*1 to fire an AP. In Fig. 4 (b) and (c), 
the potential could reach 61 because CNDs provide currents 
of enough strength. But potential of Fig. 4 (b) increase much 
faster, reach 9i by only about 4.4ms. In Fig. 4 (c), it needs 
nearly 6.7ms to fire. The reason is that for a larger M, more 
EPSPs will be generated to promote the neural firing and save 
time. 



M 


Fig. 6. fiQrx) = M ■ Qrx ■ sin(at + b), for different M and fi, T = 
100ms, 02 = 1- 

Fig. 5 presents the signal exchange process of the bit 
”1011” at the CNDs. The molecular concentration c{rk,t), the 
corresponding inspired current Ak{t) and the AP spike trains 
S{t) are respectively shown. 400ms is divided into 4 time slots, 
in which 4 bits are transmitted. In Fig. 5(a), the curve is full of 
sawteeth which are caused by noise. We could see that, during 
the time from 100 to 200 ms, low c{rk,t) means receiving little 
molecules for this bit, and the current inspired from CNDs is 
not enough to fire an AP spike. Thus this bit is ”0”. From 300 
to 400 ms, there are two AP spikes, because the molecules 
of front bit and the noise make CND receive more molecules 
than other bits, inspiring a stronger current with the amplitude 
of about 0.6 p,A/cm~^, to fire more than one spike, which 
also means one bit ” 1 ”. 

C. System Reliability 

Fig. 6 indicates the BER performance under various M and 
/i. Increasing M will propel the BER to fall first and rise then. 
The reason is two-folds. On one hand, more CNDs applied 
means more received molecules for bit ” 1 ”, correspondingly 
larger probability to fire the AP spikes in S{t). On the other 
hand, caused by channel noise and residual molecules of 
front bits, increasing CNDs quantity M will also increase the 
probability to fire the AP spikes for bit ”0”, and the spikes of 
bit ” 1 ” may also extend into bit ” 0 ” to lead the error decoding. 
Alteration of parameter p, will change the sensitivity of the 
CNDs. We could see that under larger p, fewer CNDs are 
needed to achieve the lowest BER. The result suggests that for 
different quantity of CNDs, p should be adjusted to achieve a 
lower BER. 
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Fig. 7. Pe varies with T and a, under different O 2 , M=6, ^=20nAlcm ^ 
per molecule. 

In Fig. 7, the relationship between BER, T and a under 
different 62 is revealed. We could see that, under a fixed 
setting of TX and CNDs, RX parameter 6*2 impacts BER 
directly. When 62 stays around 1, the decoding process could 
achieve minimum BER. For a longer T, the interference 
of neighbouring bit is smaller in whether diffusion-based or 
neural channel, which leads decrease of the BER. But the 
corresponding cost is a longer transmission delay. Also in 
range of 40 to 80Hz, increase the CND frequency may lead 
more AP spikes be produced and decrease the error probability, 
but higher CND working frequency requires more for nano¬ 
device manufacture. 

V. Conclusion 

In this paper, we develop a hybrid system containing 
diffusion-based and neural communication channels. We de¬ 
fine the function of the connected nano-devices, which are 
used to connect the both channels. The detail mechanism of 
exchanging information from diffusion-based to neural channel 
are proposed, and the complete communication process is 
studied in the mathematical model. The results show the 
detail process of the mechanism, indicating that it’s feasible to 
realize the communication in the hybrid system. Optimizing 
some parameters of the nano-devices could achieve lower 
BER performance for this hybrid system. Our research could 
be extended to the design of communication across various 
biological mediums in body area network. 

Although this work mainly focuses on the system model 
and information exchange mechanism of single direction, the 
future work includes the reverse direction, i.e., communication 
from neural to diffusion-based channel, as well as the analysis 
of the channel capacity, transmission delay for the hybrid 
system. 

Acknowledgment 

This work in the paper was supported by National Natural 
Science Foundation of China (NSFC-2014-61471102) and EU 
EP7 Project CLIMBER (GA-2012-318939). 


[1] I. F. Akyildiz, F. Brunetti, and C. Blazquez, “Nanonetworks: A new 
communication paradigm,” Computer Networks, vol. 52, no. 12, pp. 
2260-2279, 2008. 

[2] T. Nakano, Y. Okaie and J. Q. Liu, “Channel Model and Capacity 
Analysis of Molecular Communication with Brownian Motion,” IEEE 
Comunication Letters, vol. 16, no. 6, pp. 797-800, 2012. 

[3] D. Kilinc, O. B. Akan, “ Receiver Design for Molecular Communication,” 
IEEE Journal on Selected Areas in Communications (JSAC), 2013. 

[4] A. Guney, B. Atakan and O. B. Akan, “Mobile Ad Hoc Nanonetworks 
with Collision-Based Molecular Communication,” IEEE Transactions on 
Mobile Computing, vol. 11, no. 3, 2012. 

[5] D. Malak, and O. Akan. “A Communication Theoretical Analysis of 
Synaptic Multiple-Access Channel in Hippocampal-Cortical Neurons,” 
IEEE Transactions on communications, vol. 61, no. 6, pp. 2457-2467, 
2013. 

[6] E. Balevi and O. B. Akan, “a physical channel model of nanoscale neuro¬ 
spike communication,” IEEE Transactions on communications, vol. 61, 
pp. 1178-1187, 2013. 

[7] Q. Liu et al, “Inter-symbol Interference Analysis of Synaptic Channel in 
Molecular Communications,” IEEE International Conference on Com¬ 
munications Workshops (ICC), 2014. 

[8] F. Mesiti and 1. Balasingham, “Nanomachine-to-Neuron Communication 
Interfaces for Neuronal Stimulation at Nanoscale,” IEEE Journal on 
Selected Areas in Communications (JSAC), 2013. 

[9] B. Atakan and O. B. Akan, “Single and Multiple-Access Channel 
Capacity in Molecular Nanonetworks,” in Proc ICST/ACM Nano-Net, 
pp. 1-10, 2009. 

[10] N. Rouach, et al, “Gap junctions and connexin expression in the normal 
and pathological central nervous system,” Biology Cell, vol. 94, no. 7-8, 
pp. 457-475, 2002 

[11] J. Hjorth, K. T. Blackwell, and J. Hellgren Kotaleski, “Gap junctions 
between striatal fast-spiking interneurons regulate spiking activity and 
synchronization as a function of cortical activity,” J. Neuroscience, vol. 
29, no. 16, pp. 5276-5286, 2009. 

[12] R. D. Traub, et al, “A possible role for gap junctions in generation of 
very fast EEG oscillations preceding the onset of, and perhaps initiating, 
seizures,” Epilepsia, vol. 42, no. 2, pp. 153-170, 2001. 

[13] W. M. Kistler and W. Gerstner, “Reduction of the Hodgkin-Huxley 
Equations to a Single-Variable Threshold Model,” Neural Computation, 
V0I.9, no. 5, pp. 1015-1045, 1997. 

[14] E. P. Simoncelli, L. Paninski, J. Pillow, and O. Schwartz, “Characteri¬ 
zation of neural responses with stochastic stimuli,” in The New Cognitive 
Neurosciences, M. Gazzaniga, editor. MIT Press, 2004. 

[15] D. Malak, and O. Akan. Synaptic interference channel. IEEE Interna¬ 
tional Conference on Communications Workshops (ICC), 2013. 

[16] A. Manwani, “Information-Theoretic Analysis of Neuronal Communi¬ 
cation,” PhD thesis, California Inst of Technology, 2000. 

[17] J. S. Geelhoed, H. V. Hamelers and A. J. Stams, “Electricity-mediated 
biological hydrogen production,” Current Opinion in Microbiology, vol. 
13, issue. 3, page. 307-315, 2010. 

[18] R. Komuro and I. Balasingham, “Ejfects of Ion Channel Currents on 
Induced Action Potentials,” lEEE/EMBS Conference on Neural Engi¬ 
neering, 2013. 

[19] M. Veletic, et al, “On Spectrum Analysis for Nanomachine-to-Neuron 
Communications,” IEEE Conference on First International Black Sea, 

2013. 

[20] F. Mesiti and I. Balasingham, “Correlated Neuronal Activity in Net¬ 
works of Neurons Stimulated with Nanomachines,” ACM Nanocom, 

2014. 




